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We present measurements from BNL AGS Experiment E864 of the \H invariant multiplicity and 
of the 90% Confidence Level upper limit on the \H yield in central 11.5 A GeV/c Au + Pt collisions. 
The measurements span a rapidity range from center of mass, y C m, to j/ cm + l and a transverse 
momentum range of 0. < pr < 1-5 GeV/c. We compare these results with E864 measurements of 
stable light nuclei and particle unstable nuclei yields of the same baryon number. The implications 
""^ ' of these results for the coalescence of strange clusters are discussed. 
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Relativistic heavy ion collisions are the main experimental tool for studying the behaviour of nuclear matter un- 
der extreme energy and baryon density conditions. In addition, these collisions offer the only method to produce 
multistrange bound systems in a controlled manner, since there is copious strangeness production in these collisions. 

Hypernuclei, which are nuclei where a nucleon is replaced by a A hyperon, exist and have been studied for many 
years. More exotic forms of multi-strange nuclear systems have been hypothesized to exist. These include both 
hypernuclei which contain more than one hyperon, so called MEMOS (Metastable Exotic Hypernuclear Objects) ffl 
and strangelets [@,|[||] which are single "bags" of approximately equal numbers of strange, up and down quarks. In 
many cases the quantum numbers of the proposed MEMOS and of strangelets are the same. If so, and if strangelets 
exist, the MEMOS would decay into the more deeply bound strangelets . The production of these exotic hypernuclei 
could then be a doorway to the production of strangelets. In Experiment 864, 30 billion 10% most central, 11.5 GeV/c 
per nucleon Au on Pt or Pb collisions were sampled for the search for strangelets with A < 100 and lifetimes greater 
than 50 ns. No strangelets were observed at a level of w 10~ 8 per central collision J1m@]. 

The study of the production of the light hypernuclei 3 H and \H is essential in understanding the production 
mechanism of exotic objects such as multi-hypernuclei or MEMOS and the strangelets they might decay into if the 
latter are more stable. There are various possible production mechanisms for multihypernuclei and str ang elets such 



as the QGP distillation scenario [pPUlOjl, coalescence mechanisms |l|,[l if and thermal models jl2|,[13 14 lq| . When a 
number of nucleons coalesce to form a nucleus there is a "penalty factor" involved in the yields of the various nuclei 
for each nucleon that is added to a cluster. In the case of hypernuclei there is also an additional suppression factor 
due to the different yields of strange baryons as compared with protons and neutrons and there might also be an extra 
"strangeness penalty factor" for coalescing strange particles. The study of light nuclei in E864 jl6| is informative 
about the coalescence process of nucleons at freeze-out and the "penalty factor" involved when adding a nucleon to 
a cluster. When the invariant yields of light nuclei with A=l to A=7 are examined in a small kinematic region near 



1 



the center of mass rapidity and at low px (pt/A < 300MeV), they show an exponential dependence with baryon 
number, suggesting a penalty factor of approximately 48 for each nucleon added fig] . However in order to determine 
whether there is some extra "strangeness penalty factor" when hyperons are coalesced, the study of \H and \H 
is important. Finally the production of \H and \H in relativistic heavy ion collisions is a novel measurement and 
interesting in itself for understanding the strangeness degree of freedom in hadronic systems. In this paper we present 
measurements from BNL AGS Experiment E864 of the \H invariant multiplicity and of the 90% Confidence Level 
upper limit on the \H yield. 

II. THE E864 SPECTROMETER 

Experiment 864 is an open geometry, high data rate spectrometer designed primarily for the search of strange quark 
matter produced in relativistic Au + Pt collisions. The open geometry allows for a large region of the produced phase 
space for heavy clusters to be sampled. A beam of Au ions with momentum 11.5 GeV/c per nucleon is incident on 
a fixed Pt target. The interaction products can be identified by their charge and mass in the tracking system and 
by their energy and time of flight in the calorimeter. A detailed description of the E864 apparatus is given in Jl7| . 
Diagrams of the plan and elevation views of the apparatus are given in Figure [l]. 

The tracking system consists of two dipole analyzing magnets (Ml and M2) with vertical fields, three scintillator 
time of flight hodoscope planes (HI, H2 and H3) and two straw-tube stations (S2 and S3). The dipole magnets 
Ml and M2 can be set to different field strengths to optimize the acceptance for various particles of interest. The 
three hodoscope planes measure the time, charge and spatial position for each charged particle that passes through 
them. The straw tube planes provide improved spatial resolution for the charged tracks. The magnetic rigidity, 
momentum, and mass of the tracked particles can be determined by the position, time, and charge information from 
these detectors together with the knowledge of the magnetic field and the assumption that they come from the target. 
The hadronic calorimeter measures the energy and time of flight for all particles. The calorimeter is the primary 
detector for identifying neutral particles and can act as a powerful tool for background rejection for charged particles. 
It has excellent resolution for hadronic showers in energy (ge/E = .34/(v A E') + -035 for E in units of GeV) and time 
(at ~ 400 ps) and is described in detail in Q . There is a vacuum tank along the beam line to reduce the background 
from beam particles interacting downstream. Just upstream of the target there are beam counters and a multiplicity 
counter which are used to set a first level trigger that selects interactions according to their centrality. The calorimeter 
energy and time of flight measurements are also used to make a level-2 trigger (LET) that rejects interactions with 
no high mass particle in the spectrometer [jl9| . 

For this study the magnetic field of the spectrometer was -0.2 Tesla, which is the optimum magnetic field for the 
simultaneous acceptance of both the decay products of the hypernuclei (n~ and 3 He or 4 He ). We triggered on the 
10% most central events as defined by our multiplicity counters and an additional high mass LET trigger which was 
set for the enhancement of 3 He and 4 He ions. The LET trigger rejected interactions that did not result in any high 
mass objects in the calorimeter by a factor of approximately 60. In this way 13.5 x 10 9 10% most central collisions 
were sampled. 

III. DATA ANALYSIS 

The light hypernuclei \H and \H decay weakly into mesonic and non-mesonic channels. Their lifetimes of w 
2 x 10~ 10 sec (cr 6 cm) imply that they decay far outside the collision fireball, so we can observe them through the 
detection of their decay products. The combination of our knowledge of the magnetic field and information obtained 
by the tracking detectors of the charge, time of flight and rigidity of the decay products can provide us with accurate 
mass reconstruction of these particles |l^Jl^]. Since E864 has good particle identification for charged particles we 
have concentrated on the following mesonic channels for the \H and \H : 

\H 7T - + 3 He, branching ratio: 25% @,|l| 
\H -> 7T- + 4 He, branching ratio: 50% 

A. \H analysis 

In order to be able to reconstruct the \H signal from its decay products the "mixed event method" is used, as 
explained below. First, events that contain at least one pair of the decay products of \H are selected. The invariant 
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mass spectrum of the 3 He nuclei and the negative pions of each such event is constructed. This "same event" invariant 
mass spectrum (SE) contains a small \H signal and a large background that is due to particles that are not the decay 
products of a hypernucleus. The background shape (Bg) is obtained by constructing the invariant mass spectrum 
of uncorrelated 3 He, w~ that come from different events that contain at least one pair of the particles of interest. 
Specifically we combine the daughter particle of one type from one event with the daughter particles of the other type 
from a number of subsequent events. We have to make sure however that the mixed event spectrum does not contain 
pairs of overlapping 3 He and ir~ , which could not be found if both tracks were in the same event. This is achieved by 
requiring that the 3 He and ir~ are in different sides of the detector horizontally (magnetic bend direction), with the 
sides assigned to give optimum efficiency for simulated decays. We then simulate the shape of the hypernucleus mass 
peak by using a GEANT simulation of the decay products oi\H passing through the apparatus and reconstructing 
their invariant mass spectrum (MC). Finally we fit a linear combination of the Monte Carlo shape of the signal and 
the mixed event spectrum to the same event spectrum, 

a x (Bg) + /3x (MC) = SE. (1) 

The determination of the parameters a, (3 allow us to measure the signal either by subtracting a x (Bg) from the (SE) 
spectrum and then integrating the subtracted spectrum over the region of the expected signal or as (3 x (MC). 

As a check of this technique we looked at the p - 7r~ invariant mass spectrum from 6.5 million events to observe 
the signal from lambda decays. In Figure || we show the same event and mixed event spectra, the M.C. signal and 
the subtracted spectrum on which we overlay the M.C. signal. The agreement of the M.C. and the data is good and 
the fit suggests a signal of 7.25 a above background. We have calculated invariant multiplicities of the A in several 
rapidity and transverse momentum bins and they are found to be in good agreement with measurements from E89I 
and E877 US 

In order to reconstruct the \H signal, we have to identify its decay products. A 3 He track is defined as a charge two 
track with rapidity less than 2.7 and reconstructed mass in the range 1. — 3.4 GeV/c 2 . The efficiency of the rapidity 
cut is incorporated in the geometric acceptance for the \H and the efficiency of the mass cut is calculated to be 97% 
by using real data. A systematic error of 1% is associated with this mass cut, resulting from varying the fit parameters 
for the mass spectrum. A pion is defined as a negative particle with measured mass in the range 0. — 0.4 GeV/c 2 if 
its measured beta is less than one. Due to our finite time resolution in the hodoscopes the measured beta of a particle 
could be greater than one; any such particle with negative charge is considered to be a pion. When identifying the 
pion we avoid imposing strict beta or mass cuts which would imply strict time of flight cuts. The reason for this 
is that both of the decay particles' time of flight measurements have a common start time. Any fluctuation of the 
start time will create correlations in the masses and betas of the 3 He and ir~ which is magnified by strict cuts and 
can create structure in the same event spectrum which is not present in the mixed event spectrum. The mass cut 
used for the identification of the negative pions is loose and its efficiency is not significantly different from 100%. It 
is also required that with the apparatus divided to two horizontal sections, the 3 He be observed in one part of the 
detector and the ir~ in the other. The same event invariant mass spectrum (SE) of the 3 He nuclei and the negative 
pions that pass all the cuts is constructed. We get the background shape (Bg) by combining the daughter particle of 
one type from one event with the daughter particles of the other type from the six subsequent events. The linear fit, 
a x (Bg) + (3x (MC) = SE, yields (3 = 3.37 ± 1.67. 

The fit parameter [3 suggests a signal of 2.0 a. Figure || shows the subtracted spectrum SE — ax Bg, with the 
M.C. shape of the signal overlaid on the data. The x 2 per d.o.f. of the fit is 1.1 which implies a confidence level of 
~ 32%. If we only perform the linear fit within ten bins of where we expect our signal to be, the signal to background 
ratio does not change much and the confidence level of the fit increases to ~ 40% which gives us confidence that we 
actually have a signal. 

Assuming that the peak in the invariant mass spectrum is a signal, the invariant yield in a rapidity and pt bin 
can be calculated. Due to low statistics the yield has to be calculated in a single rapidity (y: 1.6-2.6) and transverse 
momentum (pr- 0-1.5GeV/c) bin. The reason that this specific kinematic bin was chosen was that it had reasonable 
acceptance and therefore there is not a large dependence of the yield on the assumed input model, as discussed below. 
Restricting the analysis to this kinematic region does however decrease somewhat the significance of the resulting 
signal. 

The invariant multiplicity of the J^H is calculated according to: 

Y — - x N count 

2irp~ T AyAp T N sampled x e total x g eff 

where N count is determined by the linear fit parameters as N count — (3 x (MC) , N samp i e d is the number of sampled 
events (13.5 x 10 9 ), Ay is 1, Apx 1.5 GeV/c, pr the average transverse momentum in the px bin, etotal the geometric 
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acceptance of \H including the LET, ADDMC (reconstruction efficiency), method efficiency and r/ e ff the product of 
all the other efficiencies. These efficiencies are explained in detail in the following paragraphs. 

The geometric acceptance is the fraction of \H nuclei in this kinematic bin whose decay products traverse all the 
downstream detectors. The E864 acceptance in rapidity and transverse momentum is determined by generating a 
Monte Carlo distribution oi\H particles, allowing them to decay to 3 He and ir~ and tracking the daughter particles 
using a full GEANT simulation of the experiment. The particle hit information is recorded in each of the detectors 
and then "faked" by smearing the hits according to the detector resolutions. The "faked" data is then analysed as 
the real data but with no cuts besides the fiducial ones and if both the 3 He and ir~ traverse all the detectors the \H 
rapidity and pt are reconstructed. The geometric acceptance is the ratio of these accepted \H nuclei to the generated 
ones in the same kinematic bin. The \H nuclei are generated according to a production model that is gaussian in 
rapidity with 

dN (y-y cm ) 2 ] ,„s 

— oc exp 3 

dy 2o l y 

where y cm = 1.6 is the center-of-mass rapidity and a y — 1 which is roughly estimated by folding together the rapidity 
distribution of the A f2^,^j| and the deuteron |jq| . A Boltzmann distribution is assumed for the transverse mass 
with a temperature (inverse slope) of 450 MeV. In order to study the variation of the acceptance with the assumed 
production model different widths of the gaussian distribution were used, varying between 0.7 and 1.1 and slightly 
different temperatures, including a flat distribution in pr with < pr < 3 GeV giving a systematic error of ±9%. 

The LET efficiency is the fraction of the particles of a given species that are selected by using a specific LET 
curve. This trigger efficiency can be calculated by applying the LET curve to Monte Carlo simulated showers that the 
particles of interest create in the calorimeter. For this purpose a Monte Carlo \H distribution is generated and the 
3 He decay daughters that reach the calorimeter are examined. The peak tower energy and time associated with these 
3 He nuclei are compared with the actual energy-time curve and it is determined whether the tower fired or not the 
trigger. The ratio of the number of Monte Carlo 3 He particles that fire the LET to the total number of 3 He nuclei 
that reach the calorimeter is the trigger efficiency. This method has systematic errors that depend on the simulated 
calorimeter shower response. Also time shifts in the apparatus will cause a difference between the real data and the 
theoretical curve, so a more reliable method of calculating the efficiency, especially for low efficiencies like these of 
\H , is by using the measured numbers of 3 He that did and did not fire the trigger. The number of 3 He particles 
firing the trigger is: 

N LE T = Y x N evt x R x e LET (4) 

where Y is the production rate of 3 He, N evt the number of LET triggered events, R the trigger rejection factor, €let 
the trigger efficiency. The number of 3 He particles not firing the trigger is: 

NnonLET = Y X N evt X (1 - e LET ) (5) 

so 

£let _ N LE T 

NlET + R x N non LET 

Since the LET efficiency will be different for runs with different rejection factor R, it is calculated separately for 
various runs with different rejection factors and the overall efficiency is their weighted average. 

Since it is possible for more than one track to hit the same detector element and for these tracks therefore to not 
be reconstructed, there is a track reconstruction (or ADDMC) efficiency. The ADDMC efficiency is determined by 
using Monte Carlo tracks embedded in real events. 

Finally the method efficiency is the efficiency of requiring the 3 He to be in the left part of the detector and the 
7r~ in the right. The method efficiency is the ratio of the reconstructed Monte Carlo \H in a kinematic bin after 
and before the cuts are applied. The geometric acceptance of the \H including the LET, ADDMC and method 
efficiency {e to tai = £acc x £let x ^addmc x ^method) is listed in Table |[ Other efficiencies include the detector 
efficiency, the charge two cut efficiency, the efficiency of the cuts on the % 2 distributions of the various fits applied on 
the reconstructed tracks, and the probability that the \H will interact with the target. All the efficiencies are listed 
in Table [|. Finally only \H that are reconstructed from 3 He that fired the LET can be included in our measured 
signal. This requirement reduces our counted sample of observed 3 A H by 84%. The systematic errors involved in this 
analysis are the 9% systematic error in the calculated acceptance and a 1% systematic error from the fit to the 3 He 
mass peak. The \H invariant yield, with the statistical and systematic errors combined, is listed in Table |[ 
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As an alternative method of obtaining the \H signal, we can define the ir~ as a negative particle with z component 
of the momentum being less than 3 GeV/c with all other cuts being kept the same. The effect of the strict p z cut is 
estimated by using M.C. data and is very close tol00% efficient. We then obtain the fit parameter (3 — 2.79 ± 1.69. 
The integrated signal is compatible with the results of the first method, as expected by M.C. studies. 

Because of the marginal signal, increasing the efficiency in detecting the \H is important. One way of achieving 
this is by requiring that 3 He and tt~ tracks are separated by a specific distance | Ax\ in the x direction on each of the 
detector planes instead of requiring that each decay particle is observed in a different side of the detector. However 
due to specific considerations that have to be taken into account to ensure that the same and mixed event spectra 
are treated similarly we could not achieve a much better signal to background ratio. The results from this method 
agree with the results obtained by applying strict hodoscope cuts within the statistical errors. The invariant yields 
obtained from the three different methods are within 25% of each other. The agreement of the results obtained from 
these various methods make it more probable that the peak in the invariant mass spectrum is indeed a signal and not 
just background fluctuations. 

B. \H analysis 

The first requirement for the reconstruction of the \H signal is to identify its decay products, 4 He and ir~ . A 4 He 
is defined as a charge two object with measured rapidity less than 2.7 and reconstructed mass in the range 3.2 — 6 
GeV/c 2 . The efficiency of the rapidity cut is incorporated in the geometric acceptance for the \H and the efficiency 
of the mass cut is calculated by using real data. We estimate an efficiency of 79% for the mass cut with a systematic 
error of 1% resulting from varying the parameters of the fit to the mass spectrum. A pion is defined as a negative 
charge object with measured beta greater than one or a negative particle with mass in the range 0. — 0.4 GeV/c 2 if 
its measured beta is less than one. Then it is required that the A He is observed in the left part of the detector and 
the tt~ in the right part. The same event invariant mass spectrum (SE) of the 4 He nuclei and the negative pions 
that pass all the cuts is constructed and we get the background shape (Bg) by combining the daughter particle of one 
type from one event with the daughter particles of the other type from the six subsequent events. The fit parameter 
/3 = 0.8 ±1.4 suggests that we do not have a statistically significant signal and the figure || shows the subtracted 
spectrum SE — ax Bg, with the M.C. shape of the signal overlaid on the data. 

To estimate a 90% Confidence upper limit of \H in the rapidity range 1.8 - 2.6 and pr - 1.5 GeV, the fit 
is performed by using the SE, ME and MC spectrum for pairs with a reconstructed \H rapidity and Pr in the 
above range and it yields (3 — —0.2 ± 1.3. The various efficiencies are estimated and listed in Table |J, f° r the 
rapidity range 1.8 - 2.6 and px - 1.5 GeV/c. The fit parameter f3 is negative which implies the measured signal 
((3 x M.C. — —53 ± 344) and resultant invariant yield are unphysical due to random error. The way to deal with 
that is to calculate the number Y for the invariant multiplicity according to equation |2|, though unphysical, and the 
corresponding error SY by using the fit parameters and from these the gaussian with mean Y and variance (SY) 2 is 
constructed. The physical region of the gaussian is bounded from below by 0, so the upper limit Y\ is a number such 
that the integral from to Y\ is 90% of the integral from to infinity (physical region) (^6). In this way the 90% 
Confidence Level upper limit is established as 0.42 x 10~ 4 {GeV j c)~ 2 . 

IV. RESULTS AND DISCUSSION 
A. Comparison of Hypernuclei production to non-strange nuclei 

E864 has measured the invariant yields of light stable nuclei with mass number A=l to A=7 |16[ ]. It is very useful 
to compare the yields or limits of the light hypernuclei to the yields of normal nuclei with the same mass number A 
as it would be very informative about the coalescence mechanism when strangeness is involved. In particular such a 
comparison could indicate an extra penalty factor involved in the coalescence of strangeness, if there is any. 

In the first place the \H is compared to the 3 He nucleus. Details of the measurement of 3 He can be found in [ p^[ . 
The \H yield is approximately a factor of 20 smaller than that of the 3 He in the same kinematic region. However in 
order to make a statement of whether there is an extra penalty factor when coalescing strangeness the difference in 
the production of strange and non-strange baryons should be taken into account. Therefore the relevant quantity is 

the ratio - — where Ys H , Ys He , Y\, Y p are the invariant yields of the particles in the rapidity region 1.6 - 2.6 and 

Pt region 0. - 1.5 GeV/c for \H and 3 He and 0. - 0.5 GeV/c for the p and A so as to match velocities. Furthermore 
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the yields for different species have different kinematic dependences due to collective motion. Also the efficiency of 
detecting the \H, depends strongly on the rapidity - px region. The following ratio 



R ^.(^xlxeyE.e,: (?) 

is therefore calculated, where the index i indicates various rapidity and pt bins and ti the efficiency for detecting 
the hypernucleus at each bin. Specifically, the momentum space is divided into bins of 0.5 GeV/c in px for the 3 He 
(0.f67 GeV/c for the protons and A particles to match velocities) and 0.2 units in rapidity. 

The invariant yields for protons measured by E864 |l6| are used to calculate the weighted average yields for the 
rapidity, pt bins that are used here. For the kinematic regions where there is no measurement for protons, we use 
values obtained by the following parameterization which fits our data: 

. TTtf — 0.938, r 

Y cx m t x exp[ *- — ] x [25.8 - 0.17 x (y - y cm f] (8) 

where mt is the transverse mass in units of GeV/c 2 . Similarly for the 3 He the following parameterization 

Y cx mt x exp[- mt ~ 4 Q 5 8 ° 9 ] x [8 + 17 x (y - y cm f] (9) 

is used for the bins where there was no measurement. Finally for the A the results of experiment E89f and E877 were 
used and the production model assumed p2] | was 

Y cx exp[-m t x (4.3 + 6.5 x cosh(y - y cm ) - 4.2 x (y - y cm ) 2 )]. (fO) 

The value obtained for the ratio of Equation ^ is 

Y 3H 

R =— 4 — — =0.36 ±0.26 (ff) 

which indicates that there is an extra suppression in the coalescence production of \H of about a factor of three 
compared to that of 3 He, after accounting for the different abundances of the coalescence ingredients. 

In the case of \H the 90% upper limit can be compared to the yield of 4 He. In E864 the 4 He production has been 
measured ]l6| and the following parameterization is used for the regions where there is no measurement 

m g 727 

Y cx mt x exp[ f Q — ] x [1 + 3.17 x {y - y cm f\ (12) 



The relevant ratio for the \H is 



Y X h_ 



R =---- ^ A - - — (13) 



where is the efficiency for detecting the \H and Y^ H the (unphysical) invariant yield for the \H . The ratio R and 
the relevant error of the ratio 5R are calculated and the gaussian with mean R and variance (SR) 2 is constructed. The 
physical region of the gaussian is bounded from below by 0, so the 90% confidence upper limit i?i is a number such 
that the integral from to Ri is 90% of the integral from to infinity (physical region) . Additionally the fact that 
\H has a ground state with spin J=0 and an excited state with spin J=l, which decays to the ground state, whereas 
the 4 He has only a ground state of J=0 has to be taken into account since the invariant yields of different species are 
proportional to their spin degeneracy factors (2J+1) |l6] |. So the spin-corrected 90% Confidence Level upper limit of 
the ratio is established as 

R 2 =-xRi< 0.225 (14) 
which indicates an extra suppression in the production of 4 H of at least a factor of four as compared to that of 4 He. 
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B. Implications of \H, \H results 



There is an apparent suppression in the production of the light hypernuclei as compared with the yields of non 
- strange nuclei. However before making any conclusions about an extra penalty factor involved when coalescing 
strangeness, a number of other possible reasons for this suppression should be examined. In particular both hypernuclei 
are very weakly bound and therefore could be easily destroyed in final state soft interactions, or their relatively large 
size could be a factor in their production. 



1. Finite Size Effects 

The effect of the finite size of nuclear clusters in their production has been studied by Scheibl and Heinz [£7j . Their 
coalescence model includes the dynamical expansion of the collision zone which results in correlations between the 
momenta and positions of the particles at freeze-out. The invariant spectrum of the formed clusters with mass number 
A and transverse mass M t is proportional to some effective volume V e //(A, M t ). In the saddle - point approximation 
this effective volume is proportional to the "homogeneity volume" V/ lom (m t ) of the constituent nucleons with transverse 
mass mt and transverse momentum px, where px — Pt/A with px being the transverse momentum of the nucleus. 
The Vhomimt) is accessible through HBT measurements. 

Nuclear clusters cannot be treated as pointlike particles since their rms radii are not much smaller than the 
homogeneity radii. A quantum mechanical density matrix approach allows the inclusion of finite size effects and the 
internal cluster wave function p7[ . The number of created clusters at a specific momentum is given by projecting 
the cluster's density matrix onto the constituent nucleons density matrices in the fireball at freeze-out. In the case 
of the deuteron one of the internal wave functions considered is the spherical harmonic oscillator with size parameter 
d=3.2 fm. Under various assumptions it is shown that the yield of the deuterons is identical with the classic thermal 
spectrum with only an extra quantum mechanical correction factor Cd{Rd, Pd), where Rd and Pd are the deuteron's 
momentum and center of mass space-time coordinates in the fireball rest frame. Cd provides a measure for the 
homogeneity of the nucleon phase-space around the deuteron center-of-mass coordinates. The measured deuteron 
momentum spectra do not contain information on the point of formation, so the average correction factor over the 
freeze-out hypersurface is the relevant quantity, which has a simple approximate form: 

<C d > « 1 == (15) 

where R\\(mt), Rxijnt) are the longitudinal and transverse lengths of homogeneity for the constituent nucleons. This 
means that the approximate average correction factor only depends on the ratio of the size parameter of the deuteron's 
wave function d to the radii of homogeneity of the constituent nucleons with zero transverse momentum. 

■3i 



The \H has an rms radius of « 5/m [28|] which is much bigger than the rms radius of 3 He (1.74 fm). It is therefore 
important to investigate the effect of these different sizes to the relevant yields of the nuclei. In order to make a 
rough approximation, it is assumed that the \H is a system similar to the deuteron that consists of a A and a 2 H 
and the 3 He of a proton and a 2 H. If a harmonic oscillator internal wave function is assumed for the systems with 
a size parameter for \H equal to the mean distance of the A to the cm. of the 2 H, d^ H = \J<r 2 >\d ~ 9.8/m, 

and for the 3 He equal to the mean distance of a proton to the cm. of the 2 H, d,3 He = \J <r 2 > p d ~ 2.6/m, then 
the relevant correction factors can be calculated by using equation Various AGS experiments have estimated 
the radii of homogeneity |||,||^,|3l],[32],[|| varying from 3 to 10 fm. By substituting in equation |l^ we obtain a rough 
approximation of the two correction factors. 

Ca H 

0.41 ±0.1 (16) 



C 



He 



which implies that the yield of \H as compared to that of 3 He could be a factor between two and three smaller just 
due to size effects. 

In the case of the \H, its rms radius of 2 fm [p8j is not much bigger than the A He radius (1.41 fm). Following a 
similar treatment as before, where the \H is considered to be a A bound to a triton and 4 He a proton bound to a 
triton, with size parameters, d^ h = V< r 2 >\ t ~ 3.9 fm and di He = y/ < r 2 > p t ~ 1.9 fm we have 

d H i 

« — — . (17) 
C* He 1-15 1 ' 

The size effect is not so important in the case of the \H . 
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2. Small Binding Energy 



We have previously reported that the yields of light nuclei near midrapidity and at low px are well described by the 
formula (j%) A x exp[—B/T s ] where B is the binding energy per nucleon and inverse slope T s = 5.9±l.lMeV [Q. This 
dependence cannot be explained by the coalescence or thermal models that assume a simple exponential dependence 
on the total binding energy B of the form exp[-B/T] where the temperature T of the collisions at freeze-out are in 
the order of 100-140 MeV. 

The binding energy per nucleon in the case of the \H is 0.8 MeV whereas for the 3 He it is 2.7 MeV. Therefore just 
by taking into account the exponential dependence described above the \H yield should be only 70% of that of 3 He. 

In the case of the \H which has a binding energy per nucleon of 2.5 MeV and the A He with 7 MeV/nucleon the 
effect is even bigger. The \H yield is suppressed by an extra factor of 2.2 relative to the A He yield. 

When examining the effect of the binding energy in the production of hypernuclei, it would be informative to 
compare the \H yield to that of the particle unstable nucleus A H . In E864 the yields of particle-unstable light nuclei 
have been measured |35|. The invariant yields of A H are within 50% of the 4 He yields, even though the A He has only 
a ground state of spin J=0 whereas the spin of A H is two (J=2) which means that the A H yields should be a factor 
of five greater than those of A He. If we assume feeddown from the excited states of A H (J=1,0,1) then the yields 
should be about an order of magnitude higher. Again it is not possible to explain this apparent suppression in the 
realm of the thermal and coalescence models since the mass difference between the stable and unstable nuclei is small 
compared to the collision temperatures (« lOOMeV). The spin-corrected ratio for the \H to A H is 

R = — x ^ < 1.4 (18) 

which does not indicate any extra suppression in the production of \H as compared to that of A H . Of course this 
cannot be conclusive for the existence or not of any extra suppression as we just have an upper limit for the \H. 



V. SUMMARY 



There is an apparent suppression in the production levels of \H and \H as compared with the stable light nuclei 
of the same mass number. However the size of \H and the small binding Energy of both hypernuclei seems to be an 
important factor in determining their yields. Therefore there is no clear indication of an extra penalty factor involved 
in the coalescence of strangeness and future studies are important to clarify these coalescence factors. 
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FIG. 1. The E864 spectrometer in plan and elevation views, showing the dipole magnets (Ml and M2), hodoscopes (HI, H2, 
and H3), straw tube arrays (S2 and S3) and hadronic calorimeter (CAL). The vacuum chamber is not shown in the plan view. 



10 



x 10 2 
2000 - 

1000 

x 10 : 




Same Event 



1 oooo - 



5000 - 







4000 F 



2000 



Mixed Event 




o 

4000 



2000 



M.C, 



o 



Hi 



J: 



subtracted spectrum 



- data 

- M.C. 



1-2 . . 1.4 

7T — p invariant mass 



(GeV/c) 



1.6 



11 



FIG. 2. The top panel shows the invariant mass spectrum for p - 7r~ coming from the same event. The second panel shows 
the invariant mass spectrum for p - n~ coming from different events. The following panel shows the M.C. A signal. The bottom 
panel has the subtracted invariant mass spectrum for p - tt~ . The solid histogram overlaid on the data is the M.C. A signal 
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FIG. 3. Subtracted invariant mass spectrum for 3 He - tv when strict hodoscope cuts are applied on the data. The solid 
histogram overlaid on the data is the M.C. \H signal normalized so that the peak bin matches the data. 
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FIG. 4. Subtracted invariant mass spectrum for 4 He - tv when strict hodoscope cuts are applied on the data. The solid 
histogram overlaid on the data is the M.C. %H signal 
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rapidity 
P t (GeV/c) 


1.6 - 2.6 
- 1.5 


event 


13.5 x 10 a 


count 


1220 ± 854 


ttotal = face X tLET X tADDMC X ^method 


1.96 x 1(T 4 


detector efficiency: 7/det 


0.82 2 


charge cut efficiency: r\ q 


0.84 


target absorption probability: rjtarg 


0.78 


\ 2 cut efficiency: r] x 2 


0.90 


invariant yield (GeV/c)^'' 


(5.27 ±4.04) x 10~ 4 



TABLE I. Efficiencies and invariant yields for \H in 10% most central Au-Pt collisions in units of c 2 /GeV 2 for y:1.6-2.6 and 
Pt- 0-1.5 GeV/c 



rapidity 


1.8 - 2.6 


p T (GeV/c) 


- 1.5 


event 


13.5 x 10 9 


ttotal = face X dLET X dADDMC X e met h ot J 


9.3 x 10~ 4 


detector efficiency: ?7d e t 


0.82 2 


charge cut efficiency: r?q 


0.84 


target absorption probability: rjtarg 


0.89 


\ 2 cut efficiency: rj x 2 


0.90 



TABLE II. Efficiencies for the \H in 10% most central Au-Pt collisions in units of c 2 /GeV 2 for y:l. 8-2.6 and p T : 0-1.5 
GeV/c 
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